Scaling and universality in rock fracture 
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We present a detailed statistical analysis of acoustic emission time series from laboratory rock frac- 
ture obtained from different experiments on different materials including acoustic emission controlled 
triaxial fracture and punch-through tests. In all considered cases, the waiting time distribution can 
be described by a unique scaling function indicating its universality. This scaling function is even 
indistinguishable from that for earthquakes suggesting its general validity for fracture processes 
independent of time, space and magnitude scales. 

PACS numbers: 62.20.Mk,91.30.Dk,89.75.Da,05.65.+b 



The fracture of materials is technologically of enor- 
mous interest due to its economic and human cost Q. 
Despite the large amount of experimental data and the 
considerable efforts undertaken 2] , many questions about 
fracture have not yet been answered. In particular, there 
is no comprehensive understanding of rupture phenom- 
ena but only a partial classification in restricted and 
relatively simple situations. For example, many mate- 
rial ruptures occur by a "one crack" mechanism and a 
lot of effort is being devoted to the understanding, de- 
tection and prevention of the nucleation of the crack 
0, S, H S S Q. Exceptions to the "one crack" rup- 
ture mechanism are heterogeneous materials such as fiber 
composites, rocks, concrete under compression and ma- 
terials with large distributed residual stresses. In these 
systems, failure may occur as the culmination of a pro- 
gressive damage involving complex interactions between 
multiple defects and microcracks. 

In particular, acoustic emission (AE) due to micro- 
crack growth precedes the macroscopic failure of rock 
samples under constant stress 0, ^3 ^r constant strain 
rate loading 0, This is an example of the con- 

cept of "multiple fracturing" — the coalescence of spon- 
taneously occurring microcracks leading to a catastrophic 
failure — which is thought be applicable to earthquakes 
as well 0, Q, ^J. Due to this and the similarity in 
their statistical behavior, acoustic emissions can be con- 
sidered analogous to earthquake sequences. The tempo- 
ral 0,^7], spatial 1^] and size distribution [IH of AE 
events follow a power lawj,^ust as it is commonly ob- 
served for earthquakes 19, 20]. Such power-law scaling 
can be considered indicative of self-similarity in the AE 
and earthquake source process 

The time evolution of AE and earthquake data also dis- 
play considerable differences. Laboratory rock fracture is 
dominated by a large number of foreshocks while seismic- 
ity in the Earth's crust is characterized by an abundance 
of aftershocks [2l|. Here, we show that despite this dif- 



ference the probability density function (PDF) for the 
time interval between successive events is the same in 
both cases if they are rescaled with the mean waiting 
time or equivalently with the mean rate of occurrence. In 
particular, the PDF for laboratory rock fracture neither 
depends on the specific experiment nor on the specific 
material. These observations strongly suggest a univer- 
sal character of the waiting time distribution and self- 
similarity over a wide range of activity in rock fracture. 

In Ref. (22] it was shown that the PDF of earthquake 
waiting times — without distinguishing between fore- 
shocks, main shocks or aftershocks — for different spatial 
areas, time windows and magnitude ranges can be de- 
scribed by a unique distribution if time is rescaled with 
the mean rate of seismic occurrence. It was shown in 
particular that the distribution holds from worldwide to 
local scales, for quite different tectonic environments and 
for all magnitude ranges considered. This is even true 
if the seismic rate is not stationary as during aftershock 
sequences when the rate decays according to Omori's law 
[2^. In those cases, the waiting times have to be rescaled 
by the instantaneous rate instead j22| . 

Here, we analyze in the same way the waiting times 
between AE events during periods of stationary activity 
in time series of laboratory rock fracture obtained from 
different experiments. Laboratory experiments were per- 
formed on five different materials: Flechtingen sandstone 
(Fb, porosity about 7%), Bleuerswiller sandstone (Vo, 
porosity 24%), Aue granite (Ag, porosity 1.3%), Tono 
granite (To, porosity 1.7%) and Etna Basalt (Eb, poros- 
ity 2.1%). Sandstone samples were fractured at wet con- 
dition (pore pressure 10 MPa), granite and basalt sam- 
ples at dry conditions. The formation of compaction 
bands was observed in the case of Vo sandstone and 
brittle fracture in all other cases. We investigated the 
fracture of rock samples at confining pressures in the 
range 5-100 MPa in different axial loading conditions: 
constant displacement rate of 20 /im/min (CDR), AE 
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activity feedback control of loading (AFC), as described 
in Ref. |24], and punch-through loading conditions (PT) 
[2^. The threshold level of the rate control sensor al- 
lowed varying the speed of the fault propagation by three 
orders of magnitude, i.e. from mm/s in CDR tests to 
jum/s in AE AFC tests. The newly developed data ac- 
quisition system (DaxBox made by PROKEL, Germany) 
records fully digitized waveforms (16 bit amplitude res- 
olution, 10 MHz sampling rate) in 6 Gb memory buffer, 
providing zero dead time of registration (see [2^ for a 
detailed discussion). Most importantly, the data acquisi- 
tion system allows to record AE events continuously even 
for high AE activity which is especially important for the 
analysis here. The only limitation on the shortest regis- 
tered time intervals between subsequent AE events comes 
from the finite duration of the associated signals leading 
to the possibility of overlapping AE signals, for example, 
duplets or triplets. To simplify the fully automatic pro- 
cedure of onset time picking and hypocenters location, 
only the first signal is located in a sliding window of 100 
/is. Thus, the shortest time interval in the experiments 
considered here equals 100 /is, yet the onset times of AE 
arrivals at each particular channel were determined with 
accuracy of about 0.5 /is. For each experiment, we se- 
lected one or more periods of stationary activity for our 
analysis. 

In the following, we focus on two quantities to char- 
acterize each AE event: time of occurrence and AE ad- 
justed amplitude calculated according to the procedure 
described in Ref. [22I]. While in most cases we consider 
all recorded AE events, we also study the effect of de- 
tection thresholds by only considering events with AE 
adjusted amplitude A above a certain threshold Ath- In 
both cases, the AE series is transformed into a point pro- 
cess where events occur at times U with I < i < and 
therefore, the time between successive events can be ob- 
tained as Ti = ti-^i — ti. These are the waiting times 
which are also referred to as recurrence times or interoc- 
currence times. For a given fracture experiment their 
PDF is denoted by Pe{T). 

Fig. n shows the PDF of the normalized waiting times 
6 = T / (T) for different rock fracture experiments where 
(T) is the respective mean waiting time. The excellent 
data collapse implies that P(T/(T)) does not depend on 
the particular rock fracture experiment and that we can 
write 

Pe{T) = P{T/{T)e)/{T)e. (1) 

Thus for a given fracture experiment the PDF Pe{T) 
of its AE series is determined by its mean waiting time 
{T)e = (^at — ti)/{N — 1) — or equivalently the mean 
rate Re = {T)e^ — ^nd the universal scaling function 
P{0) which can be well approximated by a Gamma dis- 
tribution 

P(i9) (X 6>-(^-^) exp i-e/B), (2) 



Table I: List of analyzed rock fracture experiments. Here, Pc 
is the confining pressure, Ath is the selected AE amplitude 
threshold, (T) is the mean waiting time and N the number of 
AE events. 
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with 7 ~ 0.8, B ^ 1.4 and the prefactor fixed by normal- 
ization 28]. Therefore, we have essentially a decreasing 
power-law with exponent about 0.2, up to the largest 
values of the argument, = T/{T) about 1, where the 
exponential factor comes into play. This is statistically 
indistinguishable (at the 2a level) from the results for 
earthquake data given in Ref. Elj, namely 7 = 0.67±0.05 
and B = 1.58 ±0.15. As shown in Fig.^ this is also con- 
firmed by the PDF of an earthquake series from Southern 
California which we have included for comparison . 

In particular, P{0) is independent of the frequency- 
amplitude distribution of the AE signal. As Fig.[21shows, 
the PDF P{A) for the AE adjusted amplitude A depends 
crucially on the particular experiment. While the overall 
structure of P{A) is rather similar for the different exper- 
iments — a sharp increase up to a maximum value fol- 
lowed by a power-law like decrease — details such as the 
location of the maximum and the slope (3 of the tail can 
be very different. For instance, the latter varies between 
2.5 and 3.7 for the considered experiments. This corre- 
sponds to a variation in the Gutenberg- Richter exponent 
b between 1.1 and 1.9 which characterizes the frequency- 
magnitude relation of earthquakes US 1^ • Thus for our 
rock fracture experiments, the variation in b is larger than 
the regional variability in the earthquake data studied 
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Figure 1: (Color online) Probability density function of nor- 
malized waiting times T/ (T) for five different rock fracture 
experiments (see Table |l] for details) and an earthquake cat- 
alog from Southern California for comparison (see text for 
details). The solid line corresponds to a fit based on Eq. (0 
with e = T/(T) giving 7 ^ 0.8 and B ^ lA. 
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Figure 2: (Color online) Probability density function of AE 
adjusted amplitude A for the same rock fracture experiments 
as in Fig. Q] 



in [221 • Yet, P{0) remains basically unchanged. More 
importantly, even considering only events above a lower 
threshold Ath-, as for the data set To5_2 in Fig.^ does not 
affect P{0). This indicates the robustness of our results. 

This robustness is further confirmed by Fig. |S1 which 
shows P{9) for a large selection of fracture experiments 
(see Tabled for details). Again, P(l9) can be well de- 
scribed by Eq. Q. The slight variation in the fitted 
values of 7 and B can be attributed mainly to statis- 
tical fluctuations and partially to measurement induced 
biases: The relatively high value of 7 for sandstone is a 
consequence of the inability to detect the shortest wait- 
ing times due to measurement restrictions absent in the 
other experiments. This absence of short waiting times 



(an order of magnitude compared with granite) signifi- 
cantly biases the estimate of 7 towards higher values. 

Fig. ini shows not only that for sandstone and differ- 
ent types of granite the influence of the specific material 
on P{9) is neglectable but also that the type of experi- 
ment (punch-through vs. constant displacement rate vs. 
activity feedback control) has no significant influence on 
P{0). Moreover, Fig. indicates that variations with Ath 
are neglectable as well. Even restricting the included AE 
events to arbitrarily selected areas within the rock sample 
did not alter P{6) (not shown). All these observations 
strongly suggest that P{0) given in Eq. Q is a universal 
result for rock fracture. It further implies that P(T) is 
self-similar over a wide range of activity rates spanning 
two orders of magnitude for the experiments considered 
here alone (see Tabled- 

Our results also indicate that the universal form of 
P{0) can be recovered for AE signals with largely vary- 
ing AE rates, as for example during foreshock sequences, 
if instantaneous rates are used. As Tabled shows, the AE 
signal of experiment Vo2 consists of at least two long sta- 
tionary regimes, Vo2_a and Vo2_b, with different (T)'s. 
Yet, the respective PDFs P{0) are indistinguishable as 
follows from Fig. |31 This implies that PiO) for the com- 
bined signal is the same as well 32j . 

While we have presented strong evidence that P{0) is 
universal for AE signals in rock fracture and earthquake 
sequences, the correlations between subsequent waiting 
times are very different. In Ref. 33], it was shown that 
the distribution of waiting times between earthquakes 
strongly depends on the previous waiting time, such that 
small and large waiting times tend to cluster in time. We 
find that this is not the case for the AE signals studied 
here. In contrast, the conditional PDF P{0\9q) is inde- 
pendent of the previous waiting time Tq with Oq = Tq/ (T) 
and, thus, P(i9|l9o) = P(<9). This might be due to the 
small number of pronounced foreshock and aftershock 
clusters of which the latter are particularly dominant in 
seismicity. 

To summarize, we have shown that the probability den- 
sity function for waiting times in laboratory rock frac- 
ture is self-similar with respect to the AE rate and can 
be described by a unique and universal scaling function 
P{0). Its particular form can be well approximated by a 
Gamma function implying a broad distribution of wait- 
ing times. This is very different from a narrow Poisson 
distribution expected for simple random processes and 
indicates the existence of a non-trivial universal mecha- 
nism in the AE generation process. The similarity with 
seismicity even suggests a connection with fracture phe- 
nomena at much larger scales and might help to under- 
stand this mechanism. 

JD would like to thank C. Goltz for stimulating dis- 
cussions. 
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Figure 3: (Color online) Probability density functions of normalized waiting times T/ (T) for fracture experiments on sandstone, 
Aue granite and Tono granite (from left to right; see Table |l] for details). The solid lines correspond to best fits according to 
Eq. (0 giving 7 = 0.92 and B = 1.2, 7 = 0.82 and B = 1.3, 7 = 0.70 and B = 1.5, respectively. 
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